Quantum phase transitions and bipartite entanglement 
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We develop a general theory of the relation between quantum phase transitions (QPTs) char- 
acterized by nonanalyticities in the energy and bipartite entanglement. We derive a functional 
relation between the matrix elements of two-particle reduced density matrices and the eigenvalues 
of general two-body Hamiltonians of d-level systems. The ground state energy eigenvalue and its 
derivatives, whose non-analyticity characterizes a QPT, are directly tied to bipartite entanglement 
measures. We show that first-order QPTs are signalled by density matrix elements themselves and 
second-order QPTs by the first derivative of density matrix elements. Our general conclusions are 
illustrated via several quantum spin models. 
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Recently, a great deal of effort has been devoted to 
the understanding of the connections between quantum 
information ^ and the theory of quantum critical phe- 
nomena 13] ■ A key novel observation is that quantum 
entanglement can play an important role in a quantum 
phase transition (QPT) ^ 1, ^ ^, i i HO, H, 

particular, for a number of spin sys- 
tems, it has been shown that QPTs are signalled by a 
critical behavior of bipartite entanglement as measured, 
for instance, in terms of the concurrence Fo'" the 

case of second-order QPTs (2QPTs), the critical point 
was found to be associated to a singularity in the deriva- 
tive of the ground state concurrence, as first illustrated, 
for the transverse field Ising chain, in Ref . |j . and gen- 
eralized in Refs. 4, 5, 6] (see Refs. HHIlIiailll for an 
analysis in terms of other entanglement measures). In the 
case of first-order QPTs (IQPTs), discontinuities in the 
gro und state concurrence were shown to detect the QPT 
|l2t IT^ The studies conducted to date are based 

on the analysis of particular many-body models. Hence 
the general connection between bipartite entanglement 
and QPTs is not yet well understood. The aim of this 
work is to discuss, in a general framework, how bipartite 
entanglement can be related to a QPT characterized by 
nonanalyticities in the energy. 

Expectation values and the reduced density matrix. — 
The most general Hamiltonian of non-identical particles, 
up to two-body interactions, reads 



{'ip\ai(3j) (7i(5j|V'), where we have used that {'ip\maif3j) 
are c- numbers and X^m 1*^) (™l ~ ^- Similarly, we can 
show that = Ti j{p^^) has a d x d matrix representa- 
tion with elements = {/3i\ p'' \ai) = (V'laj) (AlV')- 
Therefore, the energy (iplHlip) is 
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where {la.;)} is a basis for the Hilbert space, a,f3,j,5 e 
{0, 1}, and i,j enumerate N "qudits" {d- 

level systems). Let E = {ip\H\ip) be the energy 
in a non-degenerate eigenstate of the Hamilto- 
nian. The two-spin reduced density operator /5*^ is 
given by p'^ = J2m ("^IV') (V"!"^)! with m running over 
all the d^"^ orthonormal basis vectors, excluding qu- 



with U(ij) denoting a d^ x d^ matrix whose elements are 
Uap.nSj) = ^a-i^ps/^i + yafi^si whcrc TV, is the num- 
ber of qudits that qudit i interacts with, and (5^^ is the 
Kronecker symbol on qudit j. Clearly, Eq. ^ holds not 
only for the Hamiltonian operator but for any observ- 
able. Indeed, it turns out that the expectation value 
(or eigenvalue, for an eigenstate) of any two-qudit ob- 
servable in an arbitrary state is a linear function of 
the matrix elements of two-spin reduced density matri- 
ces. Moreover, it is easy to show that Eq. (0) is also 
valid for a set of qudits with distinct dimensions and 
for an arbitrary Z?-fold degenerate energy level, where 
= il/D) Y^^^i Pp , with denoting the reduced den- 
sity operator associated to the degenerate eigenstate li/^p). 
These results easily generalize to the case of a Hamil- 
tonian containing n-body terms; e.g., for a three-body 
operator O, {ip\d\^) = Eijfc Tr(0(ijA:)p*J'=), etc., for 
higher-order interactions. The above results hold for any 
value of d. Here we are especially interested in d = 2, 
i.e., the qubit case. We then use the standard basis 
{|00) , |01) , |10) , |11)} for any pair {i,j) of spins, and de- 



(O.OjIp^J \Odj), etc. 



dits i and j. 



has a d X d matrix represen- 



note p\{^{OSj\p''mj),pA 

QPT and the reduced density operator. — QPTs are 
critical changes in the properties of the ground state of 
a many-body system due to modifications in the inter- 
actions among its constituents, occurring at low tem- 
peratures T where the de Broglie thermal wavelength is 
greater than the classical correlation length of the ther- 
mal fluctuations (effectively T = 0) • Typically, such a 
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change is induced as a parameter A in the system Hamil- 
tonian H(X) is varied across a critical point Ac- Because 
they occur at T = 0, QPTs are purely driven by quantum 
fluctuations. They are associated with level crossings 
which, in many cases, lead to the presence of nonana- 
lyticities in the energy spectrum. Specifically, a IQPT is 
characterized by a finite discontinuity in the first deriva- 
tive of the ground state energy. A 2QPT (or continuous 
QPT) is similarly characterized by a finite discontinu- 
ity, or divergence, in the second derivative of the ground 
state energy, assuming the first derivative is continuous. 
These characterizations are the T = limits of the clas- 
sical definition of the corresponding phase transitions, 
given in terms of the free energy |24| . There are QPTs 
where this is not the case 0, 0| . One such example is 
the QPT in the antiferromagnetic XXZ model, where a 
critical anisotropy separates a g apless phase from a gap- 
ful phase. As shown in Ref. |23, the ground state energy 
and all of its derivatives with respect to the anisotropy 
are continuous at the critical point, despite the existence 
of the QPT. Moreover, other examples where QPTs are 
not directly related to nonanlyticities in the ground state 
energy include transitions caused by level crossings in 
the low-lying excited states [U, and those associated 
with topological order (e.g., in fractional quantum Hall 
liquids), which is not characterized by symmetry break- 
ing [23j. We shall consider in this Letter only QPTs 
characterized by nonanalytic behavior in the derivatives 
of the ground state energy. 

Assume that e^^ and V^-j^^g are smooth functions of a 
set {Xk} of couplings. If is an eigenstate of the Hamil- 
tonian then, using (9a ("0 1^/") = d\E = (-01 dxH \tp), 
we have from Eq. |(2Jl: 



(3) 



where £ = E/N. It follows immediately from Eq. Q 
that EijTr(U(ij)[9AP'^]) = 0. The origin of a IQPT 
can now be seen to be the discontinuity of one or more 
of the p^^ 's at the critical point. The second derivative, 
obtained directly from Eq. reads 



N 



aA2 



(4) 

Since U(ij) is a smooth function of {A^} and p^^ is fi- 
nite at the critical point A = Ac, it can now similarly be 
seen that the origin of the discontinuity or singularity of 
d^£{p'^^) / is due to the fact that one or more of the 
dp^^ /dX's diverge at the critical point. 

QPTs from bipartite entanglement. — In order to dis- 
cuss the role of bipartite entanglement in a QPT we 
need appropriate entanglement measures M{p^^): mono- 
tonic functions ranging from (no entanglement) to 
1 (maximal entanglement), invariant under local oper- 
ations and classical communication We consider 



two such measures: (i) concurrence C{p''^) = 

max(7j"' — ^2 ~ 73"' ~ tI"* : 0) J where the Yj ^-re the square- 
roots, in decreasing order, of the eigenvalues of the op- 
erator i?(p*'') = p^^{o-y (g) ay)p'-^*{ay (g) (Tj, ) , whcrc p*-'* 
denotes complex conjugation of p^^ in the standard ba- 
sis {|00),|01),|10),|11)}; (ii) negativity H^: J\f{p'^) ^ 
2 max(0, — mina(p^)), where are the eigenvalues of 
the partial transpose p'^^-'^^ of the density operator p'^ , 
defined as {a(3\ p-^-* |7(5) — {jf]\ p \aS). 

It is now a simple matter to connect these measures to 
the appearance of a QPT. From Eq. ^ we have E{p^^) — 
E,, Tr(U(ij)p'^) = EyTr(U^^(ij)py'^^), where the 
matrix elements of U-^-*(ij) are (a/3| U"^-* |7(S) — 
(7/3|U|a(5). Let W^^ be the unitary matrix that diag- 
onalizes p^^'^^ . Then, using Eq. we obtain 



d£{p^ 
dX 
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Theorem 1 Assume conditions (a)-(c) below are satis- 
fied. Then: a discontinuity in [discontinuity in or diver- 
gence of the first derivative of] the concurrence or nega- 
tivity is both necessary and sufficient to signal a IQPT 
[2QPT]. 

(a) The IQPT [2QPTJ is associated to a discontinuity 
in [discontinuity in or divergence of] the first [second] 
derivative of the ground state energy, which originates 
exclusively from the elements of p''^ and not, for instance, 
from the sum in Eq. ^ [Eq. ^] itself. Similarly, a dis- 
continuity in [discontinuity in or divergence of the first 
derivative of] the concurrence or negativity originates ex- 
clusively from p'^ and not from other operations such as 
max or min. 

(b) In the case of a IQPT [2QPT] the discontinuous ma- 
trix elements of p'^ present in Eq. |0j [discontinuous or 
divergent dp^^ /dX present in Eq. ^] do not either all 
accidentally vanish or cancel with other terms in the ex- 
pression for [the first derivative of] the concurrence or 
negativity. 

(c) In the case of a IQPT [2QPT] the discontinuous ma- 
trix elements of p^^ present in [discontinuous or divergent 
dp^^ / dX present in the first derivative of] the concurrence 
or negativity do not either all accidentally vanish or can- 
cel with other terms in Eq. |0j [Eq. ^]. 

Conditions (a)-(c) above are meant to exclude artifi- 
cial/accidental occurrences of non-analyticity. They are 
meant to emphasize that the entanglement-QPT connec- 
tion may directly come from the ground state reduced 
density matrix. When non-analyticities originating from 
the density operator are present in both the entanglement 
measure (or its derivatives) and the derivatives of the 
ground state energy, bipartite entanglement and QPTs 
signal each other. These observations are also the basis 
of the proof we now give. 
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Proof. IQPT: If condition (a) is satisfied then a IQPT 
must come from the discontinuity of one (or more) matrix 
elements of p'^^ , as given by Eq. |(2J). Thus, taking into 
account condition (b), the IQPT will be associated to 
a discontinuity in the concurrence or negativity, which 
is therefore a necessary condition for the IQPT. Suffi- 
ciency: (i) Concurrence - Taking into account condition 
(a), if one (or more) of the eigenvalues 7^^ of R{p^^) is 
discontinuous then one (or more) of the matrix elements 
of p^-' must be discontinuous. Assuming condition (c), a 
IQPT then follows from Eq. ||2Jl. (ii) Negativity - the 
negativity and d£{p^^) /dX are both linear in mma{p^^)- 
Therefore if the coefficient in front of miuo, (/^^ ) in Eq. @ 
does not accidentally vanish, as ensured by condition (c) , 
a discontinuous negativity signals the IQPT. 

2QPT: Considering Eq. if condition (a) is sat- 
isfied then a 2QPT must come from the discontinuity 
in or divergence of one (or more) dp^^ jdX^ since all 
the p^^ are assumed to be continuous for the case of 
a 2QPT. Thus, taking into account condition (b), the 
2QPT will be associated to a discontinuity in or di- 
vergence of the first derivative of the concurrence or 
negativity, which is therefore a necessary condition for 
the 2QPT. On the other hand, we have dxM{p'i) = 

Y.t.b=i\-(^^^(P^^)/(^Pab\^>^Pab- Therefore, taking into ac- 
count condition (a), discontinuity in or divergence of 
[dM{p^^) /dX]\^ must be caused by one or more of the 
[dp]^i,/dX]x^. Assuming condition (c), this singular be- 
havior of dp^^/dX is then a sufficient condition for a 
2QPT, which follows from Eq. Q. ■ 

Some further features following from this general anal- 
ysis are: (1) If [dM{p'^^)/dX]\^ diverges then the max- 
imal entanglement will not occur at the critical point 
Ac- (2) Concerning the behavior in the vicinity of the 
critical point: our results above show that the speed of 
divergence of both energy and the entanglement mea- 
sures is dominated by the fastest among the dp^^^^/dX 
(as illustrated in Fig. Therefore dM{p^^)/dX should 
have similar divergent properties to the second deriva- 
tive of energy. This is indeed the behavior observed for 
the transverse field Ising model in Ref. '^I. (3) Examples 
exist wherein the max/min evaluations required by the 
definition of bipartite entanglement measures generate a 
singularity related to the derivative of these measures, 
without an associated QPT condition (a) of our 

Theorem excludes such (artificial) singularities. More- 
over max/min can also eliminate singularities, a possi- 
bility which is excluded from consideration through con- 
dition (c). Next we consider examples to illustrate our 
general formalism. 

Frustrated two-leg spin-1/2 ladder. — The Hamiltonian 

— > — > ^ 

for this model is -ffiaddor = Jij S i ■ S j - hJ2 

i—l 

where Si is the spin operator vector at site i, the ex- 
change interaction along the rungs is Jij = Jr, and both 



the intra-chain nearest-neighbor and diagonal exchange 
interactions are Jjj = J. We further assume Jr > 7 J, 
with 7 « 1.401 [131. This model is exactly solvable and 
exhibits IQPTs for = Jr and = Jr+'^J ■ An anal- 
ysis of pairwise entanglement for this model can be found 
in Ref. 12] . For h < h^, and in the limit N 00, the 
ground state is a tensor product of (entangled) singlets, 
(|01) - |10))/V2, along the rungs. When hc^ < /i < hc2, 
the ground state consists of rungs which are alternately 
in singlet and (unentangled) — I triplet spin config- 
urations, |00). For h > hc2, the ground state is a tensor 
product of all rungs in the = 1 triplet state. The 
density matrix elements of the rungs are characterized 
by the following step-function discontinuities at the two 
critical points: 

n ri Vi ri f 2^ h ^ 

P22 - P33 = -P32 = -P23 = 1 0, h>h,^ 

( 0, h < hc^ , , 

= [i, h> 

where r^, with i = 1,2, denotes rungs that transition 
to the = 1 configuration at the critical point ha- 
All other density matrix elements for the rungs vanish. 
The ground state of the system is two-fold degenerate 
when hci < h < ■ The density operator for a rung is 
then represented by a statistical mixture of the broken- 
symmetry states p^^ and p'^^ , with equal probabilities. 
Indeed, for a general value of h, we can write the rung 
density matrix as p^ — {p^^ + p^^)/2. Below hc2, the 
ground state energy is given by the sum of the energies 
of each rung, due to the fact that all couplings propor- 
tional to J vanish when acting on a singlet. Using Eq. (0) 
the energy density can be then written, for h < /i^a i as 

£ = U'^iPll-P22-ph+P44 + ^iPh + P23))-HPll-Pl4)]- 

For h > hc2, contributions of the J sector must be con- 
sidered in the expression above. However the quantity 
d£/dh, which characterizes the IQPTs in this model, 
can be obtained directly from Eq. Q for any h, result- 
ing in dfiS — —\pii = —\{Pi\ + Pii)) where we have 
used that p\i^ — Q. It then follows from Eq. jnj that 
dh£ is discontinuous at both hc^ and hc2- The same 
discontinuous behavior is immediately revealed in the bi- 
partite entanglement of the spins sharing a rung. For 
these pairs a direct calculation shows that the negativ- 
ity and concurrence (which here turn out to be equal) 
read M = C = I - p\^ = I - \{p[\ + p\\), which, 
therefore, are discontinuous functions at both h = hc^ 
and h = hc2 ■ We thus find the remarkably simple result 
d£/dh — {JV — l)/2, which can also be seen as a general 
consequence of Eq. (0). This expression exemplifies how 
entanglement directly detects a IQPT. 

Permutation invariance and the transverse field Ising 
chain. — We consider now the case of Hamiltoni- 
ans whose ground states are invariant under a per- 
mutation Pik of an arbitrary pair (i, k) of spins. 
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In this case Pik = ± IV') so that (V'l = 
PkiPij PijPki- Therefore, from the general ex- 
pression for the two-spin reduced density operator 
P'^- ili^jlp'^laiPj) = X)m (7»<^jw|'0) (V'|wai/3j) = 
Y,milkSim\^p) {tl^lmakPi) = (7fci5/| p*^' |afc/3/), i.e., p'^ = 
p'^K If only a constant nearest-neighbor interaction 
is taken into account then XJ{i,i + 1) = U(j,j + 
1) = U Then, denoting Pi.i+i = Pjj+i = P 

{yi,j), we have £{p) ~ Tr(Up). As a specific exam- 
ple, consider the transverse field Ising chain with con- 
stant nearest-neighbor interactions, whose Hamiltonian 
is H ^ — J J2f=ii^'^i '^i+i + ^f)i where N is the num- 
ber of spins along the chain, af are the Pauli opera- 
tors for a spin at site i, and we use periodic bound- 
ary conditions. Setting J — 1, we obtain from Eq. Q 
that £{p) — — (i/'l cr^ — 2A(pi4 -I- P23), where the 
site-independent ground state expectation value of is 
{ip\ \ih) — pii — P44. This model presents a 2QPT at 
Ac = 1 H- This can be identified within our framework 
by taking the second derivative of £{p), yielding 
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FIG. 2: Finite size scaling of d\p22 with the number A*' of 
sites for the transverse field Ising chain. d\p22 is a function 
of N^^'^{X — Am) only, with the Ising model critical expo- 
nent p = 1, and Am being the position of the maximum of 
p22{N, A). All the data from iV = 50 to A'' = 3050 collapse 
onto a single curve. Inset: d\p22 before scaling, showing in- 
crease in singular peak sharpness with A'^ and shift of Am- 
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where we used Y.^j Tr(U(ij)^) = and Trp = 1. We 
have calculated the pij using the standard method of 
fermionization and a Bogoliubov transformation 2]. At 
the critical point Ac, Eq. (|7|) displays a divergence in the 
limit of an infinite chain. This 2QPT originates from 
the singular behavior of d\p22 and dxpa, as shown in 
Fig.n It is clear from this figure that d\p22 is dominantly 
responsible for the divergence. 
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FIG. 1: First derivative of elements of the two-spin re- 
duced density matrix for the transverse field Ising model with 
= 1000 sites. Inset: dp22/d\ and dp44/dX diverge logarith- 
mically as a function of A^ . They are fitted by x^h In N+const, 
with 0:22 = 0.135 and 0:44 = 0.024. 

Now, concerning the ground state nearest-neighbor bi- 
partite entanglement, the global 7r-rotation invariance of 
the model about the spin z-axis (.^2-symmetry) and a de- 
tailed computation of the density matrix elements leads 
to C = A/" = 2(p4i - P22)- As shown in Ref. the 



concurrence in this case is not modified by spontaneous 
symmetry breaking. In the limit N —^ 00, {dC/dX)\\^ is 
logarithmically divergent 3J. This result is here seen to 
be a direct consequence of the singular behavior of d\p22, 
just as in the second derivative of energy, since dxpn is a 
smooth function of A. Therefore d^E /dX^ and dC/dX ex- 
hibit similar critical behavior through their dependence 
upon dp22/dX, whose finite-size scaling is shown in Fig.|21 
The conclusion of Q, that the concurrence detects the 
phase transition in the Ising model, is thus simply ex- 
plained within our framework. 
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